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The z axis magnetization of a two dimensional electron gas with contact repulsive interaction and 
in presence of a Rashba potential is computed by means of quantum field theory at second order. 
A striking effect of the pure Rashba interaction is that of hindering spin alignment along the z 
direction. Evidence of transition at critical repulsive interaction coupling constant is still found. The 
degree of magnetization, however, shows a clear dependence on the spin-orbit interaction strength. 
Furthermore, the transition to magnetized state appears to be smoothed by the presence of the 
Rashba interaction. 
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I. INTRODUCTION 

The magnetization of a homogeneous gas of fcrmions 
is a long standing problem, already considered by several 
authors during the past yearsir— . It is well known, for 
instance, how, in presence of a contact repulsive interac- 
tion, both the two and three dimensional spin 1/2 fermion 
gas undergo a phase transition, magnetizing upon in- 
crease of the repulsive potential coupling constant^ii. A 
first, tentative, experimental proof of itinerant ferromag- 
netism in cold fermion atoms was recently given by Jo et 
al.^ The basic mechanism underlying the phase transi- 
tion could be well understood in terms of a competition 
between kinetic energy and repulsion of fermions with 
opposite spins. For strong repulsive interactions the po- 
tential energy reduction induced by the polarization will 
compensate the corresponding increase of kinetic energy, 
favoring the appearance of a magnetization. 
Though the physical behavior of the above mentioned 
systems appears to be well understood, the recent inter- 
est in spin orbit (SO) interactions both in semiconductor 
devices and in systems of cold atoms, leads to new ques- 
tions regarding the appearance of possible unknown spin 
effects. Among the different SO interactions present in 
semiconductors, the Rashba interactioniS certainly plays 
a fundamental role, especially due to its tunable strength. 
The form of the Rashba interaction is: 

VRashba = )^{ky(Jx — kx<7y). (1) 

where kx,y represents the momentum operator and a^^y 
are Pauli matrices. This SO interaction generates in 
asymmetric quantum wells due to the existence of an elec- 
tric field perpendicular to the plane of confinementiir— . 
The confinement introduced by the very narrow quantum 
well, in turn, allows for a two dimensional treatment of 
the electrons, which are described as effectively moving 
in the x — y plane. Moreover, an external gate poten- 
tial could be applied in order to effectively modify the 
strength of the interaction M^. 

Recently, the Rashba interaction has been reproduced in 
ultracold atoms^iii^ by means of controlled laser beams 
externally applied to the system. The realization of SO 



interactions by means of artificial gauge fields certainly 
boosted the interest on SO interactions, and a number 
of theoretical studies appeared in the last few months, 
regarding SO effects in superfluid or superconducting 
statesi^TiS^. Experimentally, a quasi two dimensional cold 
atom gas can be realized by means of counterpropagating 
laser beams along the z direction with antinodes at half 
wavelength spacing b. The ultracold atom gas is a par- 
ticularly favorable system for the study of magnetization 
properties. The absence of interfering phonons, present 
in solid state systems, and the elevated degree of con- 
trol achievable, make it a valuable tool for investigating 
delicate magnetic properties, such as combined SO and 
population imbalance effects. 

Some of the most relevant questions related to the two 
dimensional fermion gas in presence of Rashba interac- 
tion are those related to its magnetization properties. 
The interplay of Rashba interaction and magnetization 
appears indeed not yet well understood. A recent Diffu- 
sion Monte Carlo simulations^ and perturbative analyti- 
cal approaches^i for the two dimensional electron gas in 
presence of both Coulomb repulsion and Rashba inter- 
action revealed negligible two-body effects on the occu- 
pation of single particle Rashba states and no apprecia- 
ble magnetization along the z axis. In the following we 
will consider an unpolarized cold atoms assembly with 
a repulsive two body interaction. Since due to the very 
low-energy only s-wave scattering is important, the inter- 
action may be modeled by a contact interaction acting 
only between particles of opposite spins, due to Pauli 
principle. Therefore, in the action of the uniform 2-D 
system confined to the x — y plane it will be described by 
a term of the form: 

ki,k2,q 

The coupling constant g is assumed to be positive, and 
indicate the fermion fields with ±1/2 z spin compo- 
nent. Experimentally, g may be modified by exploiting 
the Feshbach resonance mechanism.— Notice that the in- 
teraction in Eq. ^ differs from the Coulomb one be- 
tween electrons in a quantum well, as it lacks the long 
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range tail and acts only between particles with opposite 
spin. 

The present paper is organized as follows: In the II sec- 
tion, the path integral computation of the free energy 
at second order in the coupling constant g will be illus- 
trated. In particular, results will be derived considering 
the possible presence of an external Zeeman potential, 
introducing a chemical potential difference, i.e. an unbal- 
ance, between the f and \. populations. Some analytical 
results will be discussed regarding the case of zero ex- 
ternal potential. In Section III numerical results will be 
shown for the magnetization of the system as a function 
of both the A and the g coupling constants. 



II. SECOND ORDER PERTURBATION 
THEORY 

In the following, a path integral derivation of the sys- 
tem free energy will be given, up to second order in the 
repulsive potential coupling constant g, while exactly in- 
cluding the Rashba effects in the independent particle 
propagator. We stress that the second order approxima- 
tion to the free energy will be obtained as an expansion in 
the two body interaction, starting from the independent 
particle solutions. This approach is to be preferred in two 
dimensions with respect to a saddle point approximation. 
In fact, in absence of Rashba interaction, the stationary 
points of the action do not correspond to its minima but 
instead they provide maxima and therefore one cannot 
proceed to evaluate the partition function by considering 
only small fluctuations around them, since large fluctu- 
ations would be favored instead. A minimum of the en- 
ergy, however, is correctly recovered by a minimization 
of the energy calculated by standard perturbation theory, 
starting from the solution of the non interacting system. 
This may be understood by expressing the energy at first 
order in g in terms of the spin f and spin \, densities: 



E/V 



-{n: 



(3) 



In the 2D gas, both the kinetic energy and the repulsive 
term show quadratic dependence As a consequence, for 
a given a,t n = -\- ri]^ it will be a parabolic function of 
the net magnetization s = — n^, with a minimum or 
a maximum at s = according to whether — g is posi- 
tive or negative, respectively. So the stationary value at 
s = is a maximum for large values of g while the energy 
is obtained when either = n or — n, i.e. when the 
system is completely polarized. The independent parti- 
cle propagator Go in presence of Rashba SO, could be 
obtained by inverting the following expression for Gq^, 
written in the f, J, spin basis: 



Gq ^(w„,k) 



(4) 



Afj. 



X{ky + ikx) 



where a;„ = ^{2n + 1) arc fermions Matsubara frequen- 
cies, ^ = is the chemical potential of the system, 
while Afi = . An external potential of the form 
—Afiaz is included, favoring the occupation of f states 
with respect to J,. The diagonalization of Gq^ leads to 
the eigenenergies 



e± = ± VAm^ + A2fc2 

corresponding to the eigenstates 



(5) 



(6) 



where the normalization constants c±^\^ are defined to be 
real and obey the equation: 



„2 

c±,k 
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(A/i ± VAm^ + A2/fc2)2 
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These states will be hereafter referred to as ± for simplic- 
ity. Knowing the Gq^ eigenstates, it is possible to write 
the transformation matrix 



c+.k- 



7^k 



c+,k 



Xiky —ikj. ) 
C- k ^ ' C_ k 



(8) 



which diagonalizcs the independent particle inverse prop- 
agator, yielding 



eo:L,('^,k) = 7^keo-l(^,k)7^[ 



(9) 



where G^^i^^g is the diagonalized independent particle in- 
verse. We stress that the inclusion of the Rashba inter- 
action in the independent particle propagator allows for 
a non perturbative treatment of the SO interaction. 
A perturbative approach is instead employed for the two 
body repulsive interaction. By using the relation 



(10) 

it is possible to employ a double Hubbard- Stratonovich 
transformation after introducing the auxiliary fields p 
and (^z, associated to density and z magnetization re- 
spectively: 

.f*<....- = /.,.e-»>-.-«-*-"-*. (U, 

The grand canonical partition function can then be ex- 
pressed as 

(12) 
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Integrating over the fermion fields and expanding the ae- 
tion up to quadratic order in the fields (p^ and p one 
obtains 

S{<i^.,p) = j drdr{ g{(fl - p") - TrlnQ^^ - 
-gTr[ga{pl - 0.(t")] + 

+ ^Tr[g^{pl - 4>,a')g^{p\ - 4>za')] + 0(5^)}. (13) 

Then, by expressing in terms of Godiag aud TZ (see 
©), it is possible to rewrite the second line of the above 
equation as: 



L • *(fc = 0) 



(14) 



where $ = (p, i/i^) and 



rT_.^( /(?.+ ) + /(C) 

(15) 

Following a similar, though slightly more complex proce- 
dure, the third line of (fT3|) is expressed as: 



with: 



(16) 



(17) 



where the trace Tr\^ is taken over two fermion momenta 
ki and k2 satisfying the relation k2 — ki = — q. M stands 
for 

Mk,,k.= E ^s,,s,.MMQk:Gl7 (18) 

Si,S2 = ± 

with similar expressions for J\f and V. Qq and Qq are 
the + and - components of Godiag and n, m and p are 
defined as follows: 



(Csi.kiAfci) 



(Ap + siy/AjJ+y^y 

(Cs2,k2Afc2)^ 



/ A m — T . „ , o . „ (Csi,kiCs2,k2) 

,S2 ,A;2 ^ '^si,S2-ki.k2 ~^ ^l,si.S2:ki .k2 ^2,si ,S2 ,^2 

Psi,S2,ki,k2 ~ ^Si,S2,fci,fc2 '^l,Si,S2,fcl,fc2 '^2,51,82,^1,^2 



(19) 



with 



C^l,Si,S2,fel,fe2 — ^('^Sl ,kl '^S2,k2 ) 

(Csi,kiCs2,k2'^)^ 



-^2, SI ,S2,fcl ,/C2 



(A^ + Siv/AAi2 + A2fc2) 

ki • k2 



(20) 



The last term on the first line of ([T3|) could also be ac- 
counted for by modifying the matrix Q into 



. 2— 2^^ ^ (21) 



The partition function (|12p can thus be written as 

= e-^«(i?etQ')"'^'e^L"2'"L (22) 

where 5*0 represents the action of the independent parti- 
cle system. The logarithm of the determinant of Q' will 
then be expanded in g and only terms up to quadratic 
order will be retained. The above expression contains 
a multiplicity of terms and certainly appears more com- 
plicated than the corresponding formula, obtained in ab- 
sence of Rashba interaction. The reason of this complica- 
tion resides in the momentum dependent transformations 
TZu , which needs to be taken into account due to the spin 
structure of the Rashba independent particle solutions. 



III. POLARIZATION 

The polarization along the z axis defined as 



M = 



(23) 



where n-f^i are the density components with spin ^ or ^ 
and n is the total density, is simply 2(0^)/n with: 



(/ D4)''Dpe-'^+-''^''(j)'') 



(/ D<f)^Dpe- 



-^In ( / D^^Dper^'+^'l''^^ 



(24) 



The introduction of the field J corresponds to modifying 
the vector L into 



L' =L 



(25) 



which then yields 
1 



[(0,l)Q'-iL + L^Q' 
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1 r fg' 
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gpV - Y^rM^ 



E/(^k)(i-24,k) + .mk)(i-2^,k) 



(26) 
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From the above formula it folfows that if A/i = then 
c±_k = l/v^ and the z polarization, is identically zero. 

We stress that this result is not sufficient in order to 
establish the absence of magnetization for any value of 
the coupling constant 5, since it is a direct consequence 
of the fact that the expectation value of (cr^) on both 
our + and — unperturbed Rashba states (see Eq. dH)) is 
zero, when = 0. Therefore even when the occupation 
numbers of these states are different, the total magne- 
tization is zero. To allow a nonzero magnetization, we 
must start from unperturbed Rashba states where (ct^) 
may be nonzero. This is easily achieved by introducing a 
variational field favoring the occupation of f over \, spin 
states, i.e. substituting A/x with + /i, in all the above 
formulas. A the end of the calculation one will get h by 
minimizing E — hMn . 

As will be discussed in the subsequent section, in gen- 
eral one will find that when Afj. = 0, the minimum energy 
is achieved for h = only when the repulsion strength g 
is below a critical value gc- No magnetization will thus 
be present when g < gc ■ 

As mentioned above, this property appears to be closely 
related to the single particle spin properties of the sys- 
tem: in fact, when /i = 0, the expectation value of <Tz 
over any of the Rashba states ([5]) (at A/i = 0) is equal to 
zero due to spin rotation around the particle wave vector 
axis2^. Remarkably, in this case the independent particle 
Rashba states contain no dependence on the parameter 
A, also implying constant (ct^) = even at A — >■ 0. 
In order to gain a more complete picture of the whole 
magnetization properties of the gas, calculations were 
also performed for {ax) and {dy). In fact, the absence 
of a z magnetization in presence of Rashba SO coupling 
does not in principle exclude the appearance of non zero 
in-plane magnetization. As an example, in the 2D fermi 
gas in absence of Rashba interaction, due to spin- rotation 
invariance, magnetization might equally occur along any 
direction if A/i = , while being necessarily oriented 
along z for A/i 7^ 0<21 

In order to describe the possible occurrence of in-plane 
magnetization, one may rewrite the contact repulsive in- 
teraction ((TU]) as^ 

^ j X! ^^^-^ \ \\ j (27) 

and consequently introduce two additional auxiliary 
fields and 0j, as done in resulting in the four 

dimensional analogues of L ([25]) and Q (fT7|) . Given the 
complexity of the calculations, these were only performed 
for the case A/t = at g < gc and will be discussed only 
qualitatively . The four dimensional analog of L will 
again show a single non-vanishing term, corresponding 
to the p field. The new matrix will instead only have non 
zero off diagonal elements corresponding to the coupling 
between (j)^ and p and between c/j^ and (f>y. As a con- 
sequence, both the expectation values of 4>x and (py are 
identically zero at A/i = 0. An analogous picture occurs 



in the Rashba interacting independent particle model, 
where both < ax > and < ay > average out to zero, due 
to the 2D rotational symmetry of the Fermi surface. 



IV. NUMERICAL RESULTS 

The numerical results presented in this section con- 
cern magnetization and quasipolarization properties of 
the system under study at T = 0. In the following we 
will express lengths in units of n~^/^,(the 2-D particle 
density), and energies in units of where m is the 
mass of the particles. Due to the elevated computational 
cost related to the treatment of some of the terms pro- 
portional to g^ in the free energy, the results presented 
in the following were obtained retaining only the terms 
of the action which are linear in g. In absence of SO 
interaction, the second order terms were shown to shift 
the transition to lower g ^ and to increase the order 
of the transition from the first to continuous. Since our 
system shows a smoothing of the transition already at 
first order in g due to the presence of SO coupling, we 
expect that second order terms will only shift gc, the crit- 
ical coupling strength for the transition, to lower values, 
without changing the overall features. 
As already outlined above, at A/i = no z magnetiza- 
tion is present in the system at repulsive couplings below 
gc- However, a phase transition, analogous to that found 
in absence of SO, is still present, causing the appearance 
of non zero magnetization at strong repulsion. Fig. [1] 
reports the magnetization as a function of the repulsion 
strength for different values of A. While no A depen- 
dence is appreciable at small g's, above the transition 
one observes a decrease of the magnetization at satu- 
ration by increasing A. Another relevant feature is the 
smoothing of the magnetization increase upon introduc- 
tion of SO, the smoothing being more relevant at higher 
A. The Rashba interaction is therefore expected to effec- 
tively frustrate M also at A/i ^ 0. This property, in fact, 
is confirmed by the results reported in Figs. 00 Above 
gc a competition will be present between the magnetiz- 
ing repulsion and the demagnetizing Rashba interaction. 
Magnetization will thus increase at fixed A by increasing 
g. When g > gc the variational parameter h introduced 
above becomes essential for the description of magneti- 
zation, since the energy minimum occurs in this case for 
/i 7^ 0, corresponding to partial spin alignment along the 
z axis. 

Non zero M is obtained, even below criticity, for Ap ^ 0, 
increasing with the A/i, as from Fig. [31 This result is 
clearly understandable, given the role of Ap in energet- 
ically favoring the occupation of = 1/2 spin states. 
Notice that the magnetization enhancement due to A/i 
decreases as the Rashba coupling constant A increases. 
The so-called " quasi-polarization" , defined as 



ntot 



(28) 
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and corresponding to the difference between the Rashba 
spin + and - density fraction at A/x = 0, shows no dif- 
ference from that of independent particles, computed by 
retaining only the kinetic energy terms in the system en- 
ergy. This also appears in good agreement with QMC 
results. For non zero A/i values, however, also the "quasi- 
polarization" appears to be affected by the presence of 
the repulsive interaction, showing dependences on g and 
A/x qualitatively similar to those of M. No direct rela- 
tion is observed between M and ^. In fact, non zero ^ is 
found in presence of zero M . At A/.J 7^ 0, however, non 
zero M and ^ are found, both increasing with g. Again, 
however, full quasipolarization is consistent with partial 
magnetization. 

V. CONCLUSIONS 

The z axis magnetization of the two dimensional 
fermion gas in presence of contact repulsion and Rashba 



SO interaction has been computed from quantum field 
theory. A general expression for the free energy at second 
order in g is obtained. As a result, we analytically found 
no magnetization at A/i = 0, at values of the repulsive 
interaction strength below a critical value gc- Above crit- 
icity magnetization is non zero and depends on both the 
SO strength and g. Numerical results were also given at 
A/i ^ 0, showing that the system develops z spin polar- 
ization by increasing A/i. Moreover, also in this case the 
repulsive interaction appears to enhance the tendency to 
develop z magnetization. The Rashba interaction acts 
by frustrating M, suggesting a possible application as an 
effective tool for controlling the system z polarization. 
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FIG. 1: (color online) Magnetization plotted in function of 
the repulsive interaction coupling constant at fixed Ajj, — 
for different values of A. Units are specified in the text. 
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FIG. 2: (color online) Magnetization plotted in function of 
the A SO coupling constant for different values of A/i at fixed 
g = 6.0 Units are specified in the text 
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FIG. 3: (color online) Magnetization plotted as a function of 
A/i at different A values and fixed g = 6.0. Units are specified 
in the text. 



